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Abstract 

Within a perturbative scalar QED model recently considered by 
Brodsky et al., we study how leading- twist Coulomb rescatterings af- 
fect the Drell-Yan cross section at small x = x tar get-, and compare to 
the case of deep inelastic scattering at small xb- We show that in 
the range where the transverse momentum transferred to the target 
is large compared to its minimal value ~ 0(x), Coulomb rescatter- 
ings affect the DIS cross section but not the Drell-Yan production 
rate. This illustrates that the leading-twist parton distribution func- 
tions become non-universal when cross sections which are differential 
in target-related particles are considered. 
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1. Introduction and summary 



Within the parton model, deep inelastic lepton-nucleon scattering structure 
functions have been shown to measure the probability to find in the tar- 
get nucleon a parton with longitudinal momentum fraction XBjorken = %b 
in the infinite momentum frame [1]. This result was obtained in a theory 
of pions and nucleons for the strong interaction. Since, the correct theory 
of the strong interaction has been established to be a gauge theory, QCD. 
According to QCD factorization theorems [2], at leading-twist the inclusive 
deep inelastic scattering (DIS) and Drell-Yan (DY) cross sections (in par- 
ticular) can be factorized and expressed as convolutions between quark and 
gluon distributions in the incoming hadron(s) and the partonic subprocess 
cross sections. The predictive power of factorization theorems arises from the 
statement that parton distribution functions are universal quantities, i.e. in- 
dependent of the collision. The universality of parton distributions appears 
to be supported by the data, at least up to some accuracy. Also, the quark 
distribution (in the nucleon N of momentum p) probed in DIS, 



where all fields are evaluated at equal light-cone time y + = and transverse 
position y± = 0j_, seems directly related to the nucleon light-cone wave- 
function in A + = gauge, supporting the probabilistic interpretation of the 
parton distribution functions (and hence of the DIS structure functions), as 
in the original parton model. 

But the expression (1) is incorrect in A + = gauge, i.e. the quark dis- 
tribution is not given by the (squared) nucleon light-cone wavefunction [3]. 
Roughly speaking, this is because in the Bjorken v — > oo limit, the eikonal 
coupling of the struck quark of momentum pi to the target color field N 1 
satisfies p\ • A oc uA + — > oo in all gauges, except A + = 0. More precisely, 
although the light-cone time y + between the absorption and emission of the 
virtual photon in the forward DIS amplitude vanishes as Coulomb in- 
teractions occurring in this short time interval actually modify the DIS cross 
section at leading-twist in all gauges, including the light-cone A + = gauge 
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[3]. Thus in a gauge theory, the simple identification between parton distribu- 
tion and parton probability (defined as the square of the nucleon light-cone 
wavefunction) does not hold. Although not excluded by this observation, 
the universality of parton distributions becomes much less intuitive. In this 
respect it was recently shown that single transverse spin asymmetries in semi- 
inclusive DIS appear at leading-twist [4], correcting previous statements [5]. 
This is due to the non-universality of spin-dependent parton distributions 
(in other words of the Sivers asymmetry [6]), which originates from the sub- 
tle behaviour of Wilson lines under time-reversal [7]. A possible correct 
expression in light-cone gauge for the gauge link entering the definition of 
spin-dependent parton distributions has recently been suggested [8]. 

In this context, it is important to reconsider the question of universality of 
spin-independent parton distributions. In the present work, I compare in a 
simple model the spin-independent quark distributions probed in DIS and in 
the Drell-Yan process at small values of x. I show that in the range of the 
exchanged transverse momentum k± responsible for leading-twist shadowing 
in DIS, the Coulomb rescattering corrections a priori modifying the DY cross 
section are in fact unitary. This is similar to what Bethe and Maximon found 
in the case of high energy bremsstrahlung and pair production [9]. Before 
the advent of QCD, it was also found that corrections to the parton model 
Drell-Yan formula are actually absent [10, 11]. In the context of gauge the- 
ories, our result is an example of the non-universality of the leading-twist 
parton distributions, which arises when considering a cross section which is 
differential in the target structure. However, according to the QCD factoriza- 
tion theorem for inclusive cross sections, we expect the k ^-integrated quark 
distributions probed in DIS and DY to be identical, even though the typical 
k± contributing in both cases is different. We check in Appendix A that this 
identity indeed occurs within a model where the scale ~ O (x) is screened by 
a finite photon mass. But the fact that this holds in general, for any target, 
is not obvious, and we think further studies are needed to settle (or disprove) 
the universality of parton distributions. 

I briefly review in section 2 the model of Brodsky et al. developed in Ref. [3] 
for DIS shadowing at small xb- We recall that this model concentrates on 
the leading-twist shadowing correction to the DIS cross section (arising from 
the aligned-jet kinematic region), which can be interpreted as part of the 
target quark distribution function probed in DIS. The typical value (k±)Dis 
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of the exchanged transverse momentum is found to be of the order of a soft 
but Xe-independent scale, (k±)ois ~ 0(m). The model is simply extended 
to DY production in section 3. Similarly to DIS, the leading-twist Coulomb 
corrections to the DY cross section arise from a kinematic region which we 
call the 'aligned-photon' region (by analogy with the aligned-jet region of 
DIS) where the longitudinal momentum fraction taken from the incoming 
projectile (anti)quark by the radiated virtual photon approaches unity. Those 
corrections are interpreted as part of the quark distribution probed in the 
DY process. We find that for Mx <C m, where M is the target mass and 
x = Xtarget <C 1, the DY cross section is unaffected by Coulomb rescattering 
at values kj_ ~ O (m), contrary to the DIS cross section. This is the main 
result of the present paper. 

This result is obtained in a scalar QED model and in the limit iC 1, which 
allows great technical simplications in the loop calculations. Since we neglect 
the scale Mx compared to kj_ from the beginning, the /cj_-integrated DY cross 
section is out of reach in the present model. Thus we cannot exclude that the 
total DY cross section receives a non-zero leading-twist shadowing correction. 
However, if this happens, the typical value of k± responsible for this effect 
must be, for x <C 1, 

(k ± ) DY ~Mi« (k ± ) DIS (2) 

This might have some implications on the properties of momentum broa- 
dening and energy loss in the Drell-Yan process. We note that the observed 
difference between the nuclear broadening of the average transverse momen- 
tum in DY production and in dijet photoproduction is not understood [12]. 
The result (2) might give some hint to this problem. 

But more importantly, it might question the universality of parton distribu- 
tions at small x, as we will discuss in section 4. In this respect, let us note 
that our result, namely the fact that Coulomb rescatterings do not modify 
the leading-twist DY Born cross section in the region of transverse momen- 
tum exchange k± ^> Mx, is similar to what was found in Ref. [13]. There it 
was shown, for transverse momenta being large compared to infrared cut-offs 
(and much smaller than the collision energy), that long-distance contribu- 
tions to the DY cross section cancel out at the two-loop order. This was 
argued to be a good indication for the validity of factorization. I stress 
that it makes on the contrary factorization much less evident, since in the 
same transverse momentum domain, Coulomb rescatterings modify the DIS 



3 



cross section, resulting in the observed nuclear shadowing of the DIS parton 
distributions. 



2. Leading- twist shadowing in DIS 

2.1 Model for the quark distribution function 

A perturbative model for leading-twist DIS shadowing has recently been stud- 
ied in [3] . Before extending this model to the DY process in the next section, 
we recall its main features. A specific contribution to ctdis is evaluated, via 
the optical theorem, from the forward DIS amplitude shown in Fig. 1. 




T(p) T(p) 

Figure 1: Forward 7*T — > 7*T amplitude in the DIS model of Ref. [3]. 

The model is perturbative and chosen to be scalar QED. One takes for the tar- 
get a scalar "quark" of mass M and momentum p, and for the light "quark" 
and "antiquark" scalars of mass m and momenta p± and p2- The couplings of 
the "gluons" of momenta hi and of the incoming virtual photon of momen- 
tum q to the scalars are denoted by g and e respectively. The forward ampli- 
tude of Fig. 1 contributes to ctdis through three different cuts between the 
Coulomb gluon exchanges. Calling A, B and C the single, double, and three- 
gluon exchange amplitudes for the process 7*(<?)T(p) — > <l(j3i)<l(p2)T(p'), the 
rescattering correction of order e 2 g 8 to the Born term / |v4| 2 reads 

Feynman diagrams contributing to A are shown in Fig. 2. 
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T(p) T(p') 

(a) (b) (c) 



Figure 2: Single gluon exchange DIS amplitude in scalar QED. 

The amplitudes B and C are obtained by adding to the Born amplitude A 
one or two gluon exchanges between the target and the light quarks. In the 
Bjorken limit 1 and at small x B , Audis receives a leading-twist contribution, 
arising from the aligned-jet configuration and presenting the features of a 
shadowing correction to the DIS Born cross section [3]. It was shown that 
the kinematic region where leading-twist shadowing appears reads: 

2v ~ V\ > P2 > k±, k i± , p 2 ±, m, M > kf, k + , p+ ~ O (Mx B ) (4) 

where when v — > oo the total momentum transfer k satisfies 

k + = Mx B + vt ( 5 ) 



The kinematic limit (4) holds in the target rest frame, where in the four- 
momentum notation k = (k + , k~ , k±) we have: 

p = (M,M,0 ± ) 

q = (-Mx B ,2u,0 ± ) 

e L = ^(1,-1, ± ) (6) 

In the case of scalar QED, the leading- twist contribution to aois arises from 
the light quarks coupling to a photon with longitudinal polarization e^. 

The scale v is the single hard scale in the problem, and the limit v — > oo 
is taken from the beginning. In the aligned-jet kinematics q~ = 1v ~ p{ , 
Coulomb rescattering corrections contribute at leading-twist to the DIS cross 

1 The Bjorken limit is defined as q~ = 2v — > oo, Q 2 = —q 2 — * oo with xb = being 
fixed. We use the light-cone variables q = q° ± q z . 
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section [14]. Compared to the scale is, the antiquark has a soft momentum p 2 
and must be considered as part of the (soft) target dynamics [2]. (At small 
xb, P2 oc 1/xb can however become large enough, so that the physics of 
destructive interferences between diffract ive amplitudes takes place, resulting 
in shadowing). In addition, the hard vertex 7*q — > q (as viewed in the infinite 
momentum frame) is taken at zeroth order in the strong coupling g. Hence 
the contribution to Aa^is arising from the domain (4) is a perturbative 
model for the scaling target quark distribution /q/T^s)- The leading- twist 
contribution to Aaois found in [3] is thus interpreted as shadowing of the 
quark distribution function in the target. 

In order to compare the quark distributions probed in DIS and DY, we will 
apply the model described above to the DY process in the next section. Let 
us repeat before the results obtained in [3] for the DIS amplitudes A, B, C 
and for the shadowing correction A<r DIS . 

2.2 DIS rescattering amplitudes 
Born amplitude 

At small Xb the Born amplitude for the DIS process is obtained in Feynman 
gauge from the dominant diagrams of Figs. 2a and 2b, and in light-cone 
A + = gauge from the diagram of Fig. 2a only. The gauge invariant result 
reads in momentum space: 



A(P2 ,P2±,k±) 



2eg 2 MQp 2 



1 



1 



(7) 



D(P2±) 



D(p 2 ± - k ± ) 



where 




2 , 2 
p ± + m n 

P 2 Mxb + rn 2 



(8) 
(9) 



In transverse coordinate space we have 



A(p 2 ,f±,R±) = 




ir± ■ P2± + iR± ■ k± 




2eg 2 MQp 2 V(m lir± )W(r ± , R ± ) 
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The functions V and W stand respectively for the incoming photon wave- 
function describing its qq content and for the qq dipole scattering amplitude: 



f d 2 p± e^- 1 '^ 1 - 
V(mr±) = J 7^y2-2 I ^ 2 = iz K ^ mr ^) 



d 2 k L 1 - 



(11) 



2 k\ 



Two-gluon exchange 

The gauge invariant expression of the one-loop DIS amplitude B correspon- 
ding to two gluon exchanges between the target and the light quarks is [3]: 



x 



B(p 2 ,P2±, k±) = -ieg 4 MQp, 
1 1 



d 2 k 



1_L 



1 



2 ' (27C) 2 k 2 ± k 2 ± 

1 



1 



D(p 2± ) D(p 2± -k 1± ) D{p 2L -k 2i _) D(p 2± -k ± ) 

— * — * — * 

where k 2 ± = k±_ — fcu_. In transverse coordinate space: 



(13) 



B(p 2 ,f ± ,R ± ) = -ieg*MQ P 2V(m\\r ± )W 2 (r ± ,R ± ) 



2! 



WA (14) 



Three-gluon exchange 

We give the expression of the three-gluon exchange amplitude C found in [3] : 

ni - - T \ 1 §-\/rr\ - f d 2 ki ± d 2 k 2 ± \ 
C{p 2 ,p 2 ±, k±) = --eg°MQp, 



2 ■ (2rc) 2 (2n) 2 k 2 1L k 2 2± k 2 



3± 



X 



+ 



1 



D(p 2 ±) D(p 2± -k 1A _) D(p 2± - k 1± - k 2± ) D(p 2± -k ± ) 



(15) 



where k 3 ± = kj_ — ki± — k 2 ±. In coordinate space: 

1 



ig2)2 W 2 A (16) 



C(p 2 ,f ± ,R ± ) = --eg 6 MQp 2 V(m liri _)W 3 (r ± ,R ± ) = 
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2.3 The Arrange in DIS 



We stress here that the amplitudes B and C are infrared finite. This is 
because the quark p 1 and antiquark p 2 form a dipole, whose scattering am- 
plitude W vanishes with the separation r± between the two quarks (see (12)). 
Thus in (13) and (15) the typical values of ki± are ~ O (k±,p 2 ±). The only 
other (soft) scale present being m\\ given in (9), the typical value of the to- 
tal exchanged transverse momentum k± contributing to the /c^-integrated 
correction Acr^is is ; 

(k±) D is ~ m\\ ~ O (m) (17) 
The rescattering correction can be obtained from Eqs. (14) and (16): 

Aa DIS ~ J d 2 f ± d 2 R ± [\B\ 2 + 2AC] = ~J d 2 f± d 2 R ± ^-A 2 W 2 (18) 

This is the leading-twist shadowing correction to the Born DIS cross section 
found in [3], interpreted as part of the (scalar) quark distribution }\/t{.xb)- 



3. Rescattering effects in Drell-Yan produc- 
tion 

3.1 Model for Drell-Yan production 

We now extend the model presented previously for DIS to the Drell-Yan 
process. This can be done by simply exchanging the virtual photon q and 
the quark p 1 . We thus describe DY production in the target rest frame where 
the incoming antiquark has a large 'minus' momentum component, p\ ~ 1v. 
As we will see the basic process for DY production in this frame corresponds 
to quark-antiquark annihilation in the infinite momentum frame. 

One diagram contributing to the DY forward amplitude is represented in 
Fig. 3. All diagrams are simply obtained by taking into account all possible 
permutations of the lower and upper vertices. 
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Figure 3: Forward amplitude of order e 2 g 8 for the Drell-Yan process. Only 
one diagram is shown. 

The Born DY cross section will get a rescattering correction: 

AaDY ~ / [ |Bmf + 2Re( ^M (19) 

where Ady, Bdy, Cdy are the amplitudes for the process q(pi)T(p) — > 
1* ( c L)Q.(P2)T {p') corresponding to one, two, and three-gluon exchange. In 
the following we will evaluate these amplitudes in the small x limit. 

In the present DY case the photon momentum q is time-like, q 2 = Q 2 > is 
the final lepton pair invariant mass squared, and the momenta are chosen as 
(<?+ > 0): " 

q = (+Mx,q~,q±) 

p = (M,M,0 ± ) (20) 

where 

* = Sr ( 21 ) 

It is easy to check that the configuration pj~ = 2is ~ q~ — > oo, which we call 
the 'aligned-photon' configuration by analogy to the DIS aligned-jet region, 
gives a leading-twist contribution to A<j DY - In the DY calculation the same 
longitudinal photon polarization vector as for DIS can be used. 

In the v — > oo limit the total momentum transfer k still satisfies (see (5)): 

k + = Mx + vt (22) 

The relevant kinematics in the target rest frame is similar to (4), 

2u = p~ > ^2 > k±, ki±, p 2 ±, k~, m, M, q ± > k? , k + = Mx+p% (23) 
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where one just added the soft q± scale. 

As in DIS, the antiquark p 2 is part of the soft target dynamics. The incoming 
"hadron" is modelled as a single antiquark, whose energy v is transferred to- 
tally to the virtual photon. Thus in the present model the colliding partons 
from the projectile and target carry respectively the momentum fractions 
x\ = 1 ~ xf and x 2 = (k + — pt)/p + = x = Q 2 / (2Mu). In the infinite mo- 
mentum frame, we recover quark-antiquark annihilation as the basic partonic 
process for DY production. 

The hard qq — > 7* vertex is still taken at zeroth order in g, thus all the 
soft dynamics should be interpreted as part of the target quark distribu- 
tion, probed at a value x of the longitudinal momentum fraction. Since the 
shadowing contribution found in DIS describes the target quark distribution 
probed at xb = x, one would naively expect, assuming parton distributions 
to be universal, to find a rescattering correction to the DY Born cross section 
originating from the domain (23) equal to that of DIS. 

As we will show, in the region (23) the rescattering corrections to a DY are 
unitary, i.e., do not modify the Born DY cross section, contrary to the DIS 
case. In this sense the effect of shifting the outgoing quark of DIS to an 
incoming antiquark in DY is drastic. 



3.2 DY rescattering amplitudes 

I now give the DY amplitudes in the small x limit. The calculation has 
been performed both in Feynman and light-cone A + = gauge, yielding 
gauge-invariant results. Since different diagrams can contribute in these two 
gauges, for simplicity the following discussion refers to the Feynman gauge 
calculation. 



Born amplitude 

The Born amplitude for the DY process is given in Feynman gauge by the 
diagrams obtained by exchanging q and p\ in Figs. 2a and 2b. The result in 
the small x limit reads: 



Ady(p 2 ,P2±,k±) 



2eg 2 MQ P2 - 
k 2 



1 



D(p 2 ±) D{p 21 _ - fcj 



(24) 



This is equal to the Born amplitude obtained for DIS (see (7)), up to an 



10 




Figure 4: Diagrams contributing dominantly to the two-gluon exchange DY 
amplitude in Feynman gauge and in the small x limit. Crossed diagrams, 
obtained by permuting the lower vertices, are included. 




Figure 5: A diagram (together with the contribution from the crossed one) 
for the two-gluon exchange DY amplitude which vanishes in the small x limit. 

irrelevant sign. This sign arises since the coupling of the photon brings a 
factor €l ■ (2pi — q) = Q for the DIS amplitude, and €l ■ (q — 2pi) for the DY 
amplitude. This is due to the fact we consider for DY an incoming antzquark 
of momentum p 1 . 

Two-gluon exchange 

In Feynman gauge, the one-loop diagrams which dominate in the small x 
limit are shown in Fig. 4. The 'crossed' diagrams, obtained by permuting 
the photon coupling vertices to the target line, are also taken into account. 
We found that the diagrams where the virtual photon emission occurs be- 
tween the two gluon exchanges are suppressed in this limit (see one example 
in Fig. 5, where the crossed diagram is also implicitly included). This sup- 
pression of radiation in DY production has been mentioned previously [15], 
but we stress here that it occurs only when the transferred momenta ki± are 
large compared to Mx, which is precisely the limit studied here (see (23)). 
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It is instructive to note the mathematical origin of this suppression, as it 
occurs in the Feynman gauge calculation. The diagram of Fig. 5 is suppressed 
because the poles in the (arbitrarily chosen) integration variable k 2 arising 
from the internal propagators p\ + k\ and P2 — k 2 lie on the same half- 
plane 2 . Note that the associated Feynman gauge DIS diagram (obtained by 
the q <-> pi exchange) is not suppressed in the small ig limit because the 
corresponding propagators read p\ — ki and p 2 — k 2 , yielding poles lying on 
different sides of the real axis. Shifting the p\ line from final (DIS) to initial 
(DY) state has non-trivial analytical consequences [13]. 

The result for the full DY one-loop amplitude is (compare to the one-loop 
DIS amplitude (13)) 



Bdy{p 2 ,P2±, k±) = -ieg 4 MQp 2 



1 



D(p 2 ±) D(p 2 ± - k ± ) 



x 



/ 



d 2 k 



1 



(27T) 2 kj ± k 2 ± 



(25) 



The infrared sensitivity of the amplitude B DY will be discussed below. 
Three-gluon exchange 

Similarly to the one-loop case, radiation within the target is suppressed in 
the region (23), where ki± ^> Mx. In Feynman gauge only two diagrams 
(including obvious permutations) contribute to the two-loop amplitude, cor- 
responding to the three exchanges occurring all before or all after the virtual 
photon emission. The result reads (compare to the DIS amplitude (15)) 



C D y{p 2 ,P2±,k ± ) 



-eg 6 MQp 2 



d 2 ki± d 2 k 



D(p 2 ±) D(p 21 _ - k ± ) 



2_L 



(2tt) 2 (2tt) 2 



P P P 

ft l_L ft 2± ft 3± 



(26) 



2 In A + = gauge the diagram of Fig. 5 can contribute, depending on the prescription 
which is used to regularize the spurious kf = pole of the gluon propagator in this gauge. 
The present discussion concerning the location of physical poles holds in any gauge, but 
in A + = gauge, a finite value for the diagram of Fig. 5 arises when one uses for instance 
the principal value prescription, because this prescription involves spurious poles on both 
sides of the real axis. 
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3.3 Absence of DY shadowing for k± ^> Mx 

We now discuss the expressions (25) and (26) for the DY loop amplitudes. 
Contrary to the case of DIS, they show an infrared sensitivity when ki± — > 0. 
This infrared singularity is absent in the cross section, the Coulomb phase 
originating from scattering between charged particles cancelling between the 
production amplitude and its conjugate. However the total DY cross section 
is out of reach within the present approximation (23). Indeed, the amplitudes 
have been evaluated with the assumption ki±,k± ^> Mx, and their precise 
infrared behaviour can thus not be inferred. However, as we will see now 
the partial contribution to the cross section originating from kj_ ^> Mx can 
be obtained. We will show that this contribution actually vanishes (at order 
eV). 

For a finite k± ^> Mx the expression (24) for the Born amplitude is valid 
and (25) and (26) can be written as 



Bdy = <7T^y / 7^ ITTT- ( 27 ) 



•t-A f d2kl± k ^ 
2 DY J (2vr) 2 k 2 1± k 2 2± 

r g 4 , r d 2 k 1± d%± kj 



One gets for the rescattering correction to the Born term: 

oc \B DY \ 2 + 2A DY C DY = \A DY \ 2 F(kl) (29) 



d 2 p 2 ±d 2 kj 

F(k 2 ± ) = ^- 4 k\ {[R 2 (k ± )] 2 - ^R 13 (k ± )] (30) 



where 



R2{ki_) 



d%_ 



k 2 1± (k±~k 1± ) 2 

12 U. . All 



Ru(k,) = ± I — d ^ i k± - (31) 

Let us stress that in the small x limit, (29) is correct for any finite k±, 
since in (29) also the momenta ki± flowing in the loops are large, ki± ^> 
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Mx. Indeed, although the individual amplitudes are infrared singular, in 
dimensional regularization one obtains the non-trivial result (see for instance 
[3] where the same expression appeared in another context): 

k ± + F(k\) oc [R 2 (k ± )} 2 - ^Ri 3 (k ± ) = (32) 

The fact that F(k\) is infrared finite shows that the typical values of ki± in 
the loop integrals of (31) are of order k±, the only scale at disposal. This 
justifies the approximation ki± ^> Mx used to evaluate the loop amplitudes. 
But since moreover F(kj_) = for any finite k±, only small k± ~ Mx — > 
may contribute to the /c^-integrated correction Actdy- 

We obtain here the main result of this paper. For a fixed k± satisfying 
k± 3> Mx, the rescattering correction (of relative order g 4 ) to the DY Born 
cross section vanishes, 

k ± » Mx =► d -^- = (33) 

This is in contrast with the DIS situation, where k± ~ mu ^> Mx contributes 
to the O (g 4 ) correction to cr DIS . These features are similar to what was found 
by Bethe and Maximon for high energy pair production and bremsstrahlung 
[9]. At momentum transfers much larger than their minimal value, Coulomb 
rescatterings modify the Born cross section for pair production, but not for 
bremsstrahlung. The absence of corrections to the parton model DY formula 
was also found in a pre-QCD context [10, 11]. 



4. Discussion 



We showed within a simple abelian model that whereas for kj_ ~ m\\ the 
DIS cross section gets a shadowing correction, Coulomb rescatterings do not 
modify the DY Born cross section at similar k±. 

It is still possible that the k±-integratedDY rescattering correction (29) could 
equal the result (18) found in DIS, in agreement with universality. But 
since the approximation (23) we used breaks down for k± ~ Mx, we cannot 
integrate (29) down to such small k±_ values and thus cannot answer this 
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question. Calculating the DY amplitudes beyond the small x limit would 
be much more involved. In particular, for k±_ ~ Mx radiation in between 
Coulomb scatterings is not suppressed. 

Since k± ^> Mx induces unitary Coulomb corrections to the Drell-Yan cross 
section, any non-vanishing contribution of order e 2 g 8 to Aa^y must arise 
from the domain k± ~ Mx, as stated in (2). The fact that different Ar- 
ranges in DIS and DY could then sum up to identical total cross sections for 
any target is not obvious. In the case of a totally screened target, with inverse 
screening length A, the values k± ~ Mx are forbidden if Mx <C A. One thus 
expects, for such values of x, Coulomb rescatterings to affect the DIS cross 
section but not the DY one (in the leading- twist regions of interest). Relying 
on Eq. (2), we thus suggest that the nucleon quark distribution functions 
probed in DIS and DY might become non-universal when M^x <C A, with 
Mtv the nucleon mass and A ~ Aqcd- We roughly estimate that the violation 
of universality sets in when the Ioffe time of the photon vjQ 2 = l/(2M^x) 
becomes larger than 1/Aqcd, i-e. when x < 0.1. 

We found instructive to supplement our scalar QED model with a mass term 
for the exchanged Coulomb photons. Calling A the photon mass, and con- 
sidering the limit A ^> Mx, the DY production amplitudes in this modified 
model are simply obtained by the replacements k\ — » k\ + A 2 in (24) and 
k 2 ± — > k 2 ± + A 2 in (25) and (26). Then (29) can be integrated over the whole 
/c_i_-range, the photon mass A acting effectively as an infrared cutoff. We 
show in Appendix A that in this specific case the integrated DY cross sec- 
tion A<j£)Y is identical (see (36)) to that of DIS given in (18). We also show 
that the typical value of k± is of order A. This illustrates that when k± can 
reach its minimal value (A in the present case), the two different Ajj_-ranges 
in DIS and DY might give equal total contributions. Let us mention that 
a similar result was found by Bethe and Maximon for pair production and 
bremsstrahlung, in the case of an unscreened target [9]. 

This somewhat academic calculation may help understanding why the uni- 
versality of the quark distribution was claimed to hold in Refs. [15, 16]. 
In these papers the DIS and DY cross sections depend on the same non- 
perturbative parameter (to be interpreted as the quark distribution, see in 
particular [15]), namely the quark pair dipole cross section in the target, ex- 
pressed in impact parameter space. It is what we find here (see the comments 
following Eq. (36)), but in the very particular case of an unscreened pointlike 



15 



target and for a finite photon mass A 3> Mx. The fact that only small k± ~ A 
contributes (see (40)) would appear difficult to infer in a coordinate space 
approach. One indeed finds that the typical value of the impact parameter in 
(36) is (R±) ~ l/ m ||- However the dominance of small k± <C m\\ for DY can 
be seen in our momentum space calculation, as expressed in (40). We explain 
in Appendix A why the relation (k±) by <S 1/ (R±) is possible. (In particular 
we do not contradict the uncertainty principle.) This point may have been 
overlooked in previous coordinate space approaches. We show in Appendix 
B that the derivation of the color dipole formulation of the Drell-Yan process 
[16, 15, 17] relies implicitly on the particular limit studied in Appendix A, 
namely A ^> Mx. Apparently no general proof, valid in the realistic limit 
A — > at fixed Mx, is known. 

The result (2) is demonstrated in the present paper by comparing leading- 
twist Coulomb rescattering corrections in DIS and DY in a model with a 
pointlike target, but which however contains the relevant features of nuclear 
shadowing [3]. Our arguments indicate that for a realistic target, leading- 
twist nuclear shadowing in DY might be reduced compared to shadowing 
in DIS. The data on DIS [18, 14] and DY [19, 20] shadowing seem to be 
reasonably consistent with the assumption that nuclear leading-twist quark 
distributions are universal, and any possible violation of universality can 
thus not be too large. But the difficulty to disentangle valence and sea quark 
shadowing, as well as quark energy loss effects [21] makes phenomenological 
analyses particularly intricate. We think that a possible violation of univer- 
sality at small x is not ruled out by the existing data. 
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Appendix A 
A particular limit: A ^> Mx 



Here we show that in the particular case where the Coulomb photons are 
given a finite mass A ^> Mx, Coulomb rescatterings affect identically the 
total DIS and DY cross sections, in agreement with universality 3 . In this 
modified scalar QED model, the DY production amplitudes (24), (25) and 
(26) are regularized in the infrared by k 2 ± — > kf ± + A 2 (and denoted by the 
subscript A), and become in transverse coordinate space: 



A^ v (p 2 ,r±,R±) 



l DY\ 



CDY(P2,r±,R±) 
G(R ± ) 



= -2eg 2 MQp 2 V(m {{ r ± ) [G(R ± ) - G(\R ± + f J) 
= -teg^MQpz V( m]l r ± ) [C(R ± ) 2 - G(\R± + f ± \ 

= +^-MQp 2 V(m\\r ± ) [G(R ± ) 3 - G(\R ± + f ± \) 



_ r d 2 k\ 
1 72tt 



(2tt) 2 k\ + A 2 2tt 



K {XR A 



(34) 



Using the above expressions one obtains 
Aa DY = 

(eg 4 MQp 2 ) 2 



d 2 f ± d 2 R ± \\Bdy\ 2 + ^dyCdy 



J d 2 r ± d 2 R ± V{m\\r 1 _) 2 [G{R ± ) - G{\R ± + r ± |)] 4 



(35) 



Assuming A <C m\\ the typical value of R± contributing to (35) is (R±) 
l/mn ~ (rj_). Using K (x) — log(l/a;) for x <C 1 we get: 



Aa DY = - l -j d 2 r ± d 2 R ± ^A 2 W 



(36) 



where A and W are given in Eqs. (10) and (12). Comparing (36) to (18) one 
sees that in this particular case (finite photon mass A ^> Mx) the leading- 

3 We stress that the target, being a scalar charged quark, is still unscreened. For a 
neutral target, the screening scale A (intrinsic to the target form factor) and the photon 
mass A would have a priori no reason to be identical. 
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twist Coulomb corrections to the total DIS and DY cross sections are iden- 
tical. Note also that similarly to [15, 16], the same quantity W 2 (the quark 
pair dipole rescattering cross section) appears in both DIS and DY results. 

We now show that whereas (k±) dis ~ m \\ contributes to the DIS cross section 
(18), (k±) DY ~ A <C m\\ contributes to (36). The equations (29), (30) and 
(31) are modified according to 



where 
R x 2 (k ± ) 

R x 13 (k ± ) 



dAa 



d 2 p 2 ±d 2 k_ 
F\k\) ee 



DY ^o,\B x 12 



'DY\ 



MA 1 2 c'A, 1.2 



l DY\ 



l F\k] 



4(2tt) 



(37) 
(38) 



d 2 k 



i± 



(*a\ + A 2 )((A; ± -A; 1± ) 2 + A 2 ) 



c? 2 A;i_l c? 2 A; 2 _l 



(*1 + A2 ) J (kj ± + X 2 )(k 2 2± + X 2 )((k ± - k 1± - k 2i _y + A 2 ) 



(39) 



For k± 3> A a difficult calculation yields: 

2 4 



fc± » A =► [R x (k ± )\ - -RUk±) ~ O ^ log 2 ^ j j (40) 

This latter equation expresses within a mass regularization scheme the result 
(32) obtained in dimensional regularization. Thus for A <C k± <C p 2 j_, the 
integrand (37) behaves as ~ \ 2 jk\ (since A^y oc k±/kj_ in this range) and 
the contribution from k±_ ^> A is suppressed in the /cj_-integrated quantity 
Aa DY - We conclude that (k±) ~ A m\\ dominates in Aa DY . In the present 
particular case (A ^ 0), this is reminiscent from the result {k±)oY ~ Mx <C 
m\\ we derived in section 3. 

The fact that (R±) ~ 1/^n and (k±) ~ A m|| does not contradict the 
uncertainty principle. One can easily see that the one-loop amplitude B^ Y 
appearing in (37) behaves as (/cj_//c 2 ) log(/c 2 /A 2 ) for A <C k± <C m\\. Hence 
/ d 2 k±\B^, Y (k±)\ 2 is dominated by the logarithmic range A <C k± <C mu. In 



18 



impact parameter space the expression / d 2 Rj_\B^ )Y (R±)\ 2 is dominated by 
the logarithmic interval 1/mii <C R± <C 1/A, as can be seen from the expres- 
sion of Bp Y given in (34), and as expected from the uncertainty principle. 
A similar conclusion is obtained for the term ~ A^yC^y in (37). However, 
performing the sum \B\ 2 + 2AC suppresses the regions k± ^> A in momentum 
space (see (40)), and R± 3> l/m\\ in coordinate space (see (35) and (36)). 
As a result (k±) ~ A and (R±) ~ V m ||- This does not contradict the un- 
certainty principle because the function \J\B\ 2 + 2AC is obviously not the 

Fourier transform of ^\B\ 2 + 2AC. We note that (k ± ) < 1/(R±) is possible 
thanks to the presence of two different scales, mu and A, and to the logarith- 
mic spread in the separate contributions from \B\ 2 and 2 AC to the DY cross 
section. This feature is absent in DIS. The DIS amplitudes are not infrared 
sensitive and thus only the scale m\\ is relevant. 



Appendix B 
Dipole formulation of the Drell-Yan process 

In this Appendix we discuss more precisely the correspondence between our 
model calculation and the dipole formulation of the Drell-Yan process [16, 15, 
17]. We first note that this formulation was originally proposed in Ref. [16], 
on the basis of a Born calculation. The Born DY diagrams obtained from 
Figs. 2a and 2b by the q <-> pi exchange were calculated in impact parameter 
space in [16]. In the present scalar QED model the Fourier transform of (24) 
reads 

A DY ( P 2, f±, R±) = -2eg 2 MQ V2 V(m {l r ± )W(f ± , R ± ) (41) 
The Born DY cross section is thus: 

°DY n ~ / d 2 r± d 2 R ± \A DY \ 2 = 4(eg 2 MQp^) 2 J d 2 f ± d 2 R ± V 2 W 2 (42) 

At the Born level, the DIS and DY cross sections are identical because the 
DIS and DY Born amplitudes given in (7) and (24) are the same (up to a 
sign). Thus the DIS quark pair dipole scattering cross section W 2 appears in 
(42). We thus recover in a simple framework the dipole formulation (obtained 
at the Born level) of the Drell-Yan process proposed in [16]. 
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In order to see what happens at higher orders, we sum over any number of 
Coulomb rescatterings. An obvious generalization of (25) and (26) yields: 



M 



DY 



Ady + Bqy + Cdy + • • • 
2ieMQp 2 ' 



D(p2±) D(p 2± - k ± ) 

1 ;„% 



d 2 ku_ 



_ l JL f 

3! J 



A(fc ± ) 



d 2 k 1± d 2 k 2 ± 



(2ir) 2 (2tt) 2 kf ± k 2 2± ki ± 



(43) 
(44) 

+ . . . 

(45) 



where k±± + . . . + k n ± = k± in the denominators k 2 ± . . . k 2 ± appearing in 
(45). As already discussed in section 3.3, the expression (45) arises from the 
assumption k±, ki± 3> Mx used to evaluate the Born and loop amplitudes, 
and contains infrared singularities when ki± — > 0. Thus integrating over 
k±, as well as Fourier transforming (44) to impact parameter space cannot 
be done, since the physical infrared regulators ~ Mx have been neglected. 
Regularizing (45) by k 2 ± — > k 2 ± + A 2 with the same parameter A in all deno- 
minators amounts to implicitly assume A ^> Mx, as was done in Appendix 
A. After this regularization is done in (45), Fourier transforming (44) gives: 



M DY (p 2 ,r ± ,R ± ) = 2ieMQp 2 V{m\\r A _) 



x 



(46) 



where G(R±) is given in (34). The equation (46) is similar to what was ob- 
tained in [15, 17]. The first and second terms arise from photon radiation 
respectively after and before the Coulomb exchanges, and involve quark im- 
pact parameters which are shifted by the amount r±. We stress here that 
these terms are infrared singular when A — > 0. Squaring (46) one gets: 



d 2 f ± d 2 R ± \Mdy\ 



(AeMQpz) 2 / d 2 f ± d 2 R ± V(m lir± y 



x sin 



-(G(R ± )-G(\R ± + f ± \)) 



(47) 



This is infrared finite when A — > 0, and equals the result obtained in DIS to 
all orders in g [3]: 

(48) 



J d 2 f ± d 2 R ± \M DY \ 2 ^ / <I-r L J-n ± | 



sm(g 2 W/2) 
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In particular the term of order g 8 in the expansion of (48) reproduces (18) 
(and (36)). The identity of the DY and DIS cross sections found here arises 
from first assuming A ^> Mx, and then taking the A — * limit. This pro- 
cedure was implicitly used in [15, 17]. In these works the step consisting in 
the replacement G(R±) - G(\R± + f±\) — >W(f ± ,R ± ) (see (47) and (48)) 

A — ^0 

is made without mentioning that G(R±) is a well-defined quantity only in 
the presence of A ^ 0. Thus the dipole formulation of the DY process 
[16, 15, 17] is in fact established beyond the Born approximation within the 
limit A ^> Mx. Independently of the question whether this formulation holds 
also in the limit A — > at fixed Mx, we stress that the very expression (46) is 
obtained from an integration over all k±s down to very small values ~ A, as 
shown in Appendix A. That such an expression can be obtained in general, 
for any realistic neutral target, has to our knowledge not been proven. 



References 

[1] S. D. Drell and T. M. Yan, Phys. Rev. Lett. 24, 181 (1970). 

S. D. Drell, D. J. Levy and T. M. Yan, Phys. Rev. D 1, 1035 (1970). 

[2] R. K. Ellis, H. Georgi, M. Machacek, H. D. Politzer and G. G. Ross, 
Nucl. Phys. B 152, 285 (1979). 

J. C. Collins and D. E. Soper, Nucl. Phys. B 194, 445 (1982). 

J. C. Collins, D. E. Soper and G. Sterman, Nucl. Phys. B 261, 104 

(1985) ; Nucl. Phys. B 308, 833 (1988); Phys. Lett. B 438, 184 (1998) 
[hep-ph/9806234], and review in Perturbative Quantum Chromodynam- 
ics (A.H. Mueller, ed., World Scientific Publ, 1989, pp. 1-91). 

G. T. Bodwin, Phys. Rev. D 31, 2616 (1985) [Erratum-ibid. D 34, 3932 

(1986) ]. 

[3] S. J. Brodsky, P. Hoyer, N. Marchal, S. Peigne and F. Sannino, hep- 
ph/0104291, to appear in Phys. Rev. D. See also S. Peigne and F. San- 
nino, Nucl. Phys. Proc. Suppl. 108, 303 (2002) [hep-ph/0112080]. 

[4] S. J. Brodsky, D. S. Hwang and I. Schmidt, Phys. Lett. B 530, 99 (2002) 
[hep-ph/0201296]. 

[5] J. C. Collins, Nucl. Phys. B 396, 161 (1993) [hep-ph/9208213]. 



21 



[6] D. W. Sivers, Phys. Rev. D 41, 83 (1990) [Annals Phys. 198, 371 (1990)]. 
D. W. Sivers, Phys. Rev. D 43, 261 (1991). 

[7] J. C. Collins, hep-ph/0204004. 

[8] X. d. Ji and F. Yuan, hep-ph/0206057. 

[9] H. A. Bethe and L. C. Maximon, Phys. Rev. 93, 768 (1954). 

[10] J. L. Cardy and G. A. Winbow, Phys. Lett. B 52, 95 (1974). 

[11] C. E. DeTar, S. D. Ellis and P. V. Landshoff, Nucl. Phys. B 87, 176 
(1975). 

[12] X. f. Guo, Phys. Rev. D 58, 114033 (1998) [hep-ph/9804234]. 

[13] W. W. Lindsay, D. A. Ross and C. T. Sachrajda, Nucl. Phys. B 214, 61 
(1983); Nucl. Phys. B 222, 189 (1983). 

[14] G. Piller and W. Weise, Phys. Rept. 330, 1 (2000) [hep-ph/9908230]. 

[15] S. J. Brodsky, A. Hebecker and E. Quack, Phys. Rev. D 55, 2584 (1997), 
hep-ph/9609384. 

[16] B. Z. Kopeliovich, Hirschegg 95, Proceedings, "Dynamical properties of 
hadrons in nuclear matter", 102-112 or [hep-ph/9609385]. 

[17] J. Raufeisen, J. C. Peng and G. C. Nayak, hep-ph/0204095. 

[18] M. Arneodo et al. [New Muon Collaboration], Nucl. Phys. B 481, 23 
(1996). 

[19] D. M. Aide et al. [FNAL E772 Collaboration], Phys. Rev. Lett. 64, 2479 
(1990). 

[20] M. A. Vasiliev et al. [FNAL E866 Collaboration], Phys. Rev. Lett. 83, 
2304 (1999) [hep-ex/9906010]. 

[21] F. Arleo, Phys. Lett. B 532, 231 (2002) [hep-ph/0201066]. 



22 



